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Abstract 

Linkod-twist maps arc area-preserving, piece-wise diffeomorphisms, 
defined on a subset of the torus. Tliey are non- uniformly liyperbolic gen- 
eralisations of the well-known Arnold Cat Map. We show that a class 
of canonical examples have polynomial decay of correlations for a-H61der 
observables, of order 1/n. 

1 Introduction 

A common method of classifying the complicated statistical properties of a dy- 
namical system is to establish its rate of decay of correlations. This is a measure 
of the rate at which the system mixes up initial conditions, independently of 
how this mixture is measured. For example, correlations for uniformly expand- 
ing maps on an interval can be easily shown to decay at exponential rate. The 
exponential nature of the decay stems, of course, from the exponential diver- 
gence of nearby initial conditions intrinsic to chaotic dynamics. Many examples 
in one dimension are now well-known, and particular interest has been shown 
to cases in which periodic boundary conditions arc replaced with an artefact 
designed to destroy the uniformity of the chaos, and hence slow the rate of 
mixing. 

Similar results in two dimensions are also established. For example, the 
Arnold Cat Map (and indeed any hyperbolic toral automorphism) can be shown 
to be exponentially mixing by appealing to the linearity of the map and using 
Fourier series [2]. This fast mixing rate has also been shown to be slowed by 
the introduction of a carefully chosen perturbation near the fixed point at the 
origin[l]. However, on the whole, interesting behaviour designed to slow down 
mixing rates tends to be restricted to behaviour at isolated points. 

In this paper we consider a linked-twist map, which could be viewed as a 
non-uniformly hyperbolic version of the miiformly hyperbolic Cat Map. It is 
Lebesgue measure-preserving, and is defined on a two dimensional manifold 
with non-trivial boundary. As such it is an instructive map, in that it reveals 



transparently both the source of its hyperbolicity, and the manner in which the 
uniformity of hyperbolicity is lost. 

The understanding of the dynamical properties of such maps was instigated 
by [10], who showed that they were almost Anosov and [3], who demonstrated 
ergodicity for a related (nonlinear) map. This was soon enhanced by [26] and 
[20], who proved mixing and the Bernoulli property respectively for families of 
linear linked twist maps. The former also treated similar examples defined on 
linked circular annuli on the plane. That these are mixing was conjectured by 
[26], the geometrical argument to demonstrate this being completed by [22]. 
(See also [21].) 

At this stage the theoretical development of such maps was left (with the 
exception of some exotic variations due to [16, 17]). However, in recent years 
[25, 23] showed that this class of maps underpins a wide variety of fluid mixing 
devices. In this context, the existence of the boundary is crucial, as it can be 
used for the first time to make rigorous statements about physically realizable 
phenomena (described in, for example [12, 11]) in practical applications to model 
the effect of hydrodynamical boundary conditions in experimental devices [24] . 
For this reason the specific details of the dynamical mechanism underlying the 
mixing properties of this particular system are likely to be of wider interest. 

In this paper however, we are concerned purely with the dynamical behaviour 
of the non-uniformly hyperbolic piece-wise diffeomorphism with boundary. We 
note that this is not the only example of a non-uniformly hyperbolic general- 
ization of the Arnold Cat map. [4] introduced another such map, also studied 
by [14] , in which non-uniformity stems from a non-monotonic twist function. In 
that case however, there exists a Markov partition which allows much immedi- 
ate analysis. In a linked twist map the dynamics are arguably more intricate, 
since a Markov partition does not exist. 

In the following we are concerned with a map on the two-dimensional torus 

= S""^ X S^. Rather than use the more standard unit interval we denote 
§^ = [0,2], with opposite ends identified (this is because we will largely be 
concerned with a subset of the torus that can now be denoted [0, 1] x [0, Ij). Let 
{x,y) S S""^ X give coordinates on T^. We define annuli 

P = X [0, 1] C and Q = [0, 1] x c T^. 

We will use the notation R = P U Q and S = P (1 Q. Define twist maps 
F : P ^ P and G : Q ^ g by 

P{x, y) = {x + 2y, y) and G{x, y) = {x,y + 2x). 

Note that F and G leave invariant the boundaries of P and Q respectively. Let 
P = id (the identity map) on R\P and G = id on R\Q so that both F and 
G are both continuous and moreover preserve the Lebesgue measure fj, on R. 
Their composition, the linked-twist map H = Go F, is illustrated in Figure 1. 
It is a Lebesgue measure-preserving piece-wise diffeomorphism of R into itself. 

[26] showed that H has the iC-property and the result of [6] shows that in fact 
it is Bernoulli. It is non-uniformly hyperbolic, the proof of non-zero Lyapunov 
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(a) P, Q and S (b) Image under F (c) Image under H 

Figure 1: Illustration of the toral linked- twist map H : R ^ R. The white region 
shows T^\i?. In (a) the region S is coloured mid-grey, with the remainder of 
P shown in dark-grey and the remainder of Q shown in light-grey. Figure (b) 
shows the image of these sets under the twist F whilst (c) shows the image 
under H = GoF. 



exponents /z-a.e. following from an invariant cone field; for details see [23]. No 
results concerning the rate of mixing for H are known to us. 

For n e N and for any pair of bounded, measurable functions (p,tp : R ^ R 
('observables') define the correlation function 

Cn{^,i;,H,n)= [ {^oH^)i;dfi- [ (/jdM / Vd/i. (1) 
Jr Jr Jr 

It is well-known that {H, ^) is mixing if and only if C„ — )• for any such pair of 

observables. 

The rate of decay of correlations for H refers to the order of this convergence 
for sufficiently regular (f and tp- Let Ha denote the space of real- valued, a-H61der 

functions on R. These arc the functions (p : R ^M. ioi which there are positive 
constants a and C so that for all z, z' d R sufficiently close 

\^{z)-^{z')\^Cd{z,zT, 

where d(-, •) denotes distance on R. As is common we make the further assump- 
tion Jj^ipd/j, = 0, which simplifies (1) at no expense of generality. 
The main result of our paper is the following. 

Theorem 1.1. Ifip,tpG Ha then \Cn{ip,ip, H, fj,)\ = 0{l/n). 

It is important to remark now that, although both statement and proof of 

Theorem 1.1 make explicit use of the particular annuli P and Q defined above, 
this restriction is little more than a notational convenience. All of our results 
hold, with only superficial alterations, in the general case P = §^ x [poiPi]) 
Q = [qo,Qi] X S^, for any choice of pq ^ Pi,9o 7^ 91 € §"'^7 with F and G 
appropriately re-defined also. We remark on this further at the end of Section 6. 

We note, given that linked twist maps can be used as a model for a wide 
variety of mixing devices [23], that Theorem 1.1 gives a practical bound on 
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mixing rates in such applications. Moreover, it can be shown that in general, the 
polynomial rate given is indeed attained by typical observables. This argument 
and its relevance in applications is discussed in detail in [24]. 

Out paper is organised as follows. Our proof of Theorem 1.1 uses certain 
recent results from the dynamical billiards literature and we give a synopsis 
of these in Section 2. The results essentially reduce the problem to a detailed 
analysis of an induced map given by first returns to S* C R. This is carried 
out in Sections 3 through 5; in particular, in Section 3 we study the partition 
of S induced by the return map, in Section 4 we show that the return map is 
Bernoulli and in Section 5 wc show that a technical condition regarding local 
expansion factors (to be defined in Section 2) is satisfied. In Section 6 we bring 
these results together to conclude the proof of Theorem 1.1. Finally in Section 7 
we collect a few thoughts regarding potential extensions and generalisations of 
our result. 

2 Decay of correlations in hyperbolic systems 

We describe some recent results concerning the decay of correlations in systems 
with some hyperbolicity, the foundations of which are to be found in two seminal 
papers of Young [27, 28] . 

Let X be a Riemannian manifold, possibly with boundary, and let T : X — >■ 
X be a hyperbolic map preserving an ergodic SRB measure u. Let A c X, of 
positive i^-measurc, have hyperbolic product structure, i.e. A is the intersection 
of a family of stable manifolds with a family of unstable manifolds. 

For X e X the first return time Rtn(x;T,A) = min{n > 1 : T"(x) e A} 
denotes the first iterate of x to enter, or return to, the sot A. Ergodicity ensures 
such a value exists almost everywhere. Of the successive returns of x, the 
first to satisfy an additional, technical condition on the length of local invariant 
manifolds (we omit the details, for which see the original papers) will be denoted 
by Rtn* (x; T, A) and called the first good return. 

Theorem 2.1 ([28]). If there exists a > such that 

lv{x e X : Rtn*(x;T,A) > n} = ©(n""), 

I then for a-Holder observables f, g we have 

lC„(/,5,T,t/)|=0(n-°). 

In practice constructing A and establishing (2.1) can be prohibitively dif- 
ficult. For the Arnold Cat Map, the procedm-e is described explicitly in [7], 
but this is a particularly straightforward construction, relying on the uniform 
hyperbolicity and linearity of the map. A more tractable method is to first find 

a set F C AT where hyperbolicity is 'strong', so that we can choose K gY G X , 
and so that the induced map Ty -.Y defined by first returns satisfies 

lv{x e X : Rtn*(a;; Ty, A) > n] = 0{9'^), 
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1 for some 9 € (0, 1). This can be achieved, without needing to exphcitly con- 
struct A, by estabUshing a few conditions first given by Chernov [5] and later 
improved upon by Chernov and Zhang [8]. These conditions are reproduced at 
the end of the present section. 

Finally (2.1) can be established from (2) by a method essentially owing to 
Markarian [15]. The method, developed a little in [8], introduces a redundant 
logarithmic factor to the decay rate, but this problem is resolved by a general 
scheme of Chernov and Zhang [9] . 

We now list the conditions given in [8] that collectively establish (2). 

Smoothness 

X is an open domain in a smooth (C°°) two-dimensional compact Rieman- 
nian manifold. The possibility of points at which T is undefined, discontinuous 
and/or non-differentiable is admitted; in this case such points are contained 
within a closed set D of zero Lebesgue measure. We refer to D as the singular- 
ity set. We denote by Dm = Ui!lo^ T~^{D) the singularity set for T™ and by 
D_rn = \jT=Q T'{D) the singularity set for T"™. 

Hyperbolicity 

There are two families of cones C"(a;) and C"'(x) in the tangent space T^X for 
X Cz X. These families are continuous on X and the angle between complemen- 
tary cones is bounded away from zero. They are invariant in the sense that 
DTi^ix)) C C"(r(a;)) and DTiCix)) D C"*(T(a;)) whenever DT exists, and 
they are expanded in the sense that 

\\DTv\\ ^ X\\v\\ for all v G C"{x)\\DT-^v\\ ^ X\\v\\ for all v G C'ix), 

where A > 1 is a constant and || • |1 the Euclidean norm. For m > 0, all tangent 
vectors to lie in stable (C*) cones and all tangent vectors to D_m lie in 
unstable (C") cones. 

If v' is a T-invariant probability measure then i^'-a.e. x G X has one positive 
and one negative Lyapunov exponent as well as one stable and one unstable 
manifold. We denote these W^{x) and W""{x) respectively. 

SRB measure 

T : X ^ X preserves a mixing measure f whose conditional distributions on 
unstable manifolds are absolutely continuous, i.e. u is an SRB measure. 

Distortion bounds 

Let X{x) denote the factor of expansion on unstable manifold a.t x E X. If 
X, y belong to the same unstable manifold and if T" is defined and smooth 
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on then 

where d(-,-) denotes distance on X and ^ : — t- M+ is some function, inde- 
pendent of the choice of W, so that ^(f) as f ^ 0. 

Bounded curvature 

The curvature of unstable manifolds is uniformly bounded by a constant B ^ 0. 

Absolute continuity 

If Wi , W2 are small, close unstable manifolds then the holonomy map h : Wi 

W2, defined (where applicable) by sliding along stable manifolds, is absolutely 
continuous with respect to the Lebesgue measures (induced by the Euclidean 
metric) on Wi and W2. Moreover the Jacobian is bounded, i.e. 

for some C > 1. Here W[ C Wi denotes those points at which h is defined. 



Structure of the singularity set 

We say that C X is an admissible curve in the unstable cone field, or more 
concisely an unstable curve, if all tangent vectors to W are in unstable cones. For 
any admissible curve W the set W OD is at most countable and has at most K 
accumulation points on W, K being a constant. Moreover if {Xn}nm dW r\D 
is a sequence converging to an accumulation point Xoo then 

d(a;„, Xoo) ^ const ■ n~'^ 

for some constant d > 0. 



One-step growth of unstable manifolds 

Let ly be a local unstable manifold, denote by Wi the connected components 
of W\D and let Aj = min{A(x) : x € Wi}, which is the minimal local expansion 
factor of T on Wj. We have 




where \W\ denotes the length of unstable manifold W and the supremum is 
taken over all unstable manifolds. The condition describes strong expansion 
along unstable manifolds. If T is not sufficiently expansive it is enough that 
T™ satisfies the condition for some m € N, with D, W and A appropriately 
redefined. 
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(b) Singularity set for Gs Singularity set for Hs, denoted cr 

Figure 2: Singularity sets for the return maps (a) Fg, (b) Gs and (c) Hs- 
Details of the constuction of singularity sets are given in appendix A. The 
dotted lines in figure (c) show the unstable manifolds of points p and q, and the 
components of the sets Si and S2 are labelled. The sets Sn for large n can be 
seen accumulating in the top-left and bottom-right corners of S (points q and p 
respectively), with some components of 6*3 , S'4 , 6*5 , S'g labelled. Note that these 
regions in the singularity set for Hs in which the Sn accumulate on p and q 
correspond to the relevant regions in the singularity sets for Fs and Gs, and so 
their structure can be easily found explicitly. 



This completes the list of conditions. 

3 The natural partition of the induced map 

To prove Theorem 1.1 we show that (2.1) is satisfied, with {H, R, ji) taking the 
place of (T, X, J/). To that end we take Y = S and establish (2) by verifying 
the conditions listed in Section 2. This occupies the present section and the 
following two. In Section 6 we establish (2.1) as described, however we do not 
need to appeal to [9] in order to avoid redundant factors, rather we can employ 
an instructive result, Lemma 3.2, of the present section. 

We begin the analysis of Hs by considering the natural partition it induces 
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on S. Let Fs : S ^ S, the return map with respect to the twist F, be given by 

Fs{z) = F"{z), where n = Rtii{z; F, S). 
Similarly define Gs : S ^ S. It is easily checked that 

Hs = Gso Fs. 

Those points 2: S 5 for which Rtn(2:; F, S) is 'large' are confined to neighbour- 
hoods of the corners p = (1, 0), q = (0, 1) e S. The same is true for Rtn(2; G, S) 
and Rtn(z; H, S). Singularity sets for Fs, Gs and Hs arc shown in Figures 2(a), 
2(b) and 2(c) respectively. Dotted lines represent local stable manifolds of p and 
q, which lie within the set of points which return to S under a single iterate 
of H. The components of the set S2 are also labelled. The structure of these 
singularity sets are given in appendix A, although the majority of our arguments 
do not require the precise geometrical details. 

We denote by a the singularity set for Hs and by cr" the singularity set for 
Hg, n € Z. The set a partitions S into countably many open sets on which Hs 
is a linear map characterised by a hyperbolic matrix 

n..-no..n.s-{l ?)'(; J / = ( ^ 

Here j = m,n{z; F, S) and k = Rtn{Fs{z); G, S). We remark that Rtn{z; H, S) = 
j + k-1. 
Let 

Sn = {zeS: Rtn{z; H, S) = n} . 

Clearly S\ U^i has zero /U-measure. Let ns denote the restriction of /x to 
S. 

Proposition 3.1. ns is an ergodic, invariant measure for Hs- 

The result is entirely standard and we omit a proof. What is not immediately 
clear is that Hs is in fact Bernoulli; we prove this in Section 4. 

For z G S let {u, v) = {dx, dy) give coordinates in the tangent space T^S" = 
M?. In T^S we define a pair of cones 

C+{z) = {{u,v) ■.u = ovv/u^ 1}, C-{z) = {{u,v) : u = or -u/v^ 1}, 

called the unstable and stable cones at z respectively. These are illustrated in 

Figure 3. The cones arc independent of the underlying point z and the angle 
between them is bounded away from zero. Define the unstable and stable cone- 
fields: 

C± = U C^iz). 

zes 

Our next result says that is invariant under and uniformly expanded by the 
derivative DHg^. Moreover contains all tangent vectors to the singularity 
set for H^^. Let || • || : — ^ [0, 00) denote the standard Euclidean norm. 



8 




C+(z) = {(u,v):u = OoTvlu>l} C-{z) = {{u,v) : v = oi -u/v> 1} 

Figure 3: The unstable and stable cones C'^{z) and C~{z). 

Lemma 3.1 (Hyperbolicity). IfwG is a tangent vector to z G S\a^^ then 

{DH^^),w C C^{H^\z))\\{DH^^),w\\ > V5\\w\\. 

Conversely, if w is a tangent vector to z G a^^ and is tangent to a singularity 
line through that point, then w G . 

Proof. For the first statement, we deal only with the case w G C+, the other 

being entirely similar. Fix z G S\a, let w = (u, v)'^ G C~^{z) be a tangent vector 
at z and le w' = {u',v')^ = DHs{u,v)^ be a corresponding tangent vector at 
Hs{z). Using (2) we have 

v' 2ku + 4:jkv 

- > Z^T^ = 2k>l, 

i.e. w' G C+{Hs{z)). Moreover 

+ w'2 > (l + 4fc2) + v"^ {Af + {Ajk + 1)^) > 5 (w^ + v^) , 

showing that > \/5||w'||- 

For the second statement, we deal only with the case z G a'^, the other 
being entirely similar. The boundary of S consists of horizontal and vertical 
lines, so if {u, is tangent to the boundary then either u = or u = 0. 
The singularity set consists of the i/g-images of the boundary, therefore 
its tangents consist of the Dils-images of those tangents vectors just described. 
(2) gives 

DHs{Q,vf = {2jv, {Ajk + l)vfDHs{u,Of = {u,2kuf, 

and it is easily shown that these vectors are in C+ as required. □ 

Our next result concerns the 'itinerary' of z G S with respect to Hs- We 
show that if n is large and z G Sn then some number, depending only on n, of 
the immediate pre- images and images of z must be in Si (recall Figure 2(c)). 
In effect long returns are isolated. This feature of the dynamics turns out to be 
crucial in establishing the polynomial decay rate. 
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Lemma 3.2 (Isolation of large return times). There are constants K,k > so 
that if z £ Sn and n > Ke^^ then H^giz) € 5*1 for each 1 ^ |i| ^ N . 

Proof. Roughly speaking, if n is large then Sn is close to p or g and Hs{Sn) is 
close to or s, respectively. The (exponential) rate at which successive images 
move away is bounded and so some number of iterates remain in ^i. The 
behaviour of Hg^ is similar. 

We deal rigorously with one case. Let 5^ C Sn be the connected component 
close to p and having a; = 1 as a boundary. Let S[ C 5i be the connected 
component adjacent to 0. By considering the map H and the partition of S in 
Figure 2(c) we observe that 

(i) there exists A > 1 such that \i z e S[ then d{{),Hs{z)) < Ad(0,z), 

(ii) there exists c > such that if z e 5*,'^ then d(0, Hs{z)) ^ c/n, 

(iii) there exists C > such that if d(z, 0) < C then z e S'J . 

Now suppose that z G S'^ and n > ■ It follows from, in turn, (i), (ii) 
and the assumption on n that 

d{0, Ws^\z)) ^ X'd{0, Hs{z)) ^ ^ CA'-^, 
and finally from (iii) that H^^^[z) G S[ for each O^i^N. □ 

4 The Bernoulli property for the induced map 

In this section we prove the following: 
Theorem 4.1. Hs is Bernoulli. 

Fundamental results concerning the ergodic properties of non-uniformly hy- 
perbolic systems were established by Pesin [19] and extended to a class of smooth 
maps with singularities by Katok et. al. [13]. We describe some results of the 
latter, restricting ourselves to the two-dimensional case. 

Let T be a map defined on an open subset X of a compact Riemannian 
manifold and preserving a measure v. T is smooth except possibly for a set of 
singularities (points of discontinuity or non-differentiability) contained within 
a union D of smooth, compact submanifolds of positive codimensions. The 
'heaviness' of D is restricted thus: there are positive constants o, Ci, and for 
any e > 

Iv {Be{D))^Cie\ 

1 Bg{D) is the ^-neighbourhood of D using the Riemannian metric. This is 
commonly refered to as the condition (KSl). A further condition (KS2) requires 
an upper bound on the growth of second derivatives in the vicinity of D, however 
it requires a little notation to give a precise formulation and will be trivially 



10 



J(t) 

e-neighbourhood — ■ 

. 2^.- 

> 2^/^ + 2e 

(a) 

Figure 4: Illustration of the e-neighbourhood for the singularity set . 

satisfied by piecewise-Unear systems such as ours, so we do not state it. The 
following condition of Oseledec [18] is also required: 

I [ log+ \\DT\\du < 00, 
Jx 

1 where log^(a;) := max{log(a;), 0} for x G X. We call it condition (OS). 

If T : X X as above satisfies (KSl), (KS2) and (OS), then Lyapunov 
exponents 

X±(x,t)= hm -log||l?r"(x)t|| 

n— >±oo n 

exist for i^-a.e. x Cz X and each tangent vector t at x. If there is one positive 
and one negative Lyapunov exponent at x then x has a local unstable manifold 
7"(a;) and a local stable manifold Y{x), and these arc absolutely continuous. 
X has an ergodic partition, meaning that X ~ \SiL\^i where T\x^ is crgodic 
for each i and, moreover, each Xi = [S^^\ Xij where T"(*)|jc.^ is Bernoulli for 
each j. 

All of these conclusions hold equally for T™, m e N. 

Lemma 4.1. Hs : S ^ S satisfies (4), i-e. there are a > 0, Ci > 0, and for 

any e > 

^^s (Be(c7)) < 

Proof. Singularity line-segments accumulate in the four 'groups' shown in Fig- 
ure 2(c). We consider one such group, shown in Figmc 4(a), and establish an 
appropriate bound. The lemma follows easily. Excluding L, index these line- 
segments cr„, n e N, in order of decreasing length. Notice that^ length((T„) ~ 
1/n; in particular, the total length is unbounded. 

Fix a small e > 0. We construct the e-neighbourhood in two stages. Let 
J(e) be a rectangle with sides parallel to coordinate directions; of height -^/e 
and width and with lower-right vertex at (1, 0), coinciding with that of S. 

^Here we use ~ to indicate the common asymptotic notation given by / ~ 5 if //g 1 
for functions / and g. 
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Sec Figure 4(a). Now let N be the smallest integer to exceed l/2y^, then J(e) 
contains cr„ for every n ^ N and 

US {B,{J{s))) = + 2e) (2\/i + 2e) = 2£ + Se^/^ ^ 4^2 ^ ^je. 

There remain N — 1 line-segments to consider. An e-neighbourhood of one 
such segment is shown in Figure 4(b). The total measure is at most 

N-l Ll/2\/^J /I \ 1 



n=l n=l 



The asymptotic notation describes the limit e — )• (equivalently A'' — )• 00) and 
wo have used the fact that X]!!Li 1/"- ^ In A^. Finally we observe that as .t ^ 00 
the polynomial x^, b > grows more quickly than Inx. Substituting x = 1/e, 
this says that as £ — >■ the polynomial e°, a < 1 grows more quickly than £ In |. 
Thus the lemma holds for any < a < 1 and an appropriate Ci. 

□ 

Lemma 4.2. Hs : S ^ S satisfies (4), i-e. 

log+ ||£>/fs||d/is < 00. 



Proof. \\DHs{z)\\ is given by the largest eigenvalue of a hyperbolic matrix as in 
(2). If z e Sn and n is large then either Hs{z) = Goi^"(0) or Hs{z) = G''oF{z) 
(this is an easy consequence of the dynamical features described in Lemma 3.2). 
In either case the eigenvalue in question is 1 + 2n + -\/4n(nr+T) ~ 4n. Thus 



log+ \\DHs{z)\\d^ls ^ const ^log+(4n)/x(5„). 

Sn (see Figure 4) is approximately rectangular, having width ~ 1/n and height 
- l/{n + 1) - 1/n = e» (l/n^) so that Ms(S'„) = O (l/n^). Thus the sum 

converges. □ 

For yus-a.e. z G S and for each non-zero w £ T^S, Lemmas 4.1 and 4.2 
establish the existence of Lyapunov exponents 

X±{z,w)= lim -log||DiJ^(^)«;|| 

ra— >±(X) n 

for Hs- Their existence for i?™, m € N, is an easy consequence. 

Lemma 4.3. Lyapunov exponents for Hs : S ^ S are non-zero. 

Proof. It is a standard result that for a /xg-typical z G S there can be at most two 
distinct Lyapunov exponents. For such a ^; let w e C^{z) C T^S. Lemma 3.1 
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shows that DH^{z)w e C"(^;) and that \\DH^{z)w\\ > 5"/2||u;||, for each n e N. 
Thus ^ ^ 

X+iz,w)> lim - log5"/^||w|| = - logS > 0, 

n— >oo n 2 

showing that z has a positive Lyapunov exponent. Similarly, taking w' € C'^{z) 
leads to the conclusion that x_(2:, w') < and thus that z also has a negative 
Lyapunov exponent. □ 

Wc now consider global properties of stable and unstable manifolds to prove 
the main theorem of this section. For clarity, we define 7" (2) and 7" (2;) to be 
local unstable and stable manifolds, respectively, of z; by definition these are 
connected. To prove theorem 4.1 wc will confirm the intersection of forward 
iterates of 7" (2;) with backward iterates of Y{z). In section 5 we will study 
global unstable and stable manifolds, defined by W^{z) = Un>o^s'T''(^) ^^'^ 
W^{z) = Un>o^s"T*(^) respectively. These are only piecewise connected as 
they are cut under iteration as described in the following proof. 

We note that the corresponding manifolds for manifolds for H can be com- 
puted explicitly, as described in [26], as straight lines with gradient given by a 
continued fraction whose entries are given by successive return times to S. A 
similar approach could be taken for Hs, but we only require to observe that 
gradients of such straight lines are constrained to lie within the cones C~^{z) 
and C~{z). 

Proof of Theorem 4-1- It is enough ([13] or [20]) that for yu-a.e. z,z' £ S and all 
m, n G N large enough 

lH^r{z)nHs"'Y{z')j^t 

1 This property is clearly satisfied by H. As discussed in [20, 26], images of local 
unstable and stable manifolds, under iteration of H, diverge exponentially in 
length, remain connected, and are mutually transversal. As soon as they span 
S in the vertical and horizontal directions respectively, the analogue of (4) is 
guaranteed. 

The difference in the argument for Hs is that although images of local in- 
variant manifolds diverge exponentially in length and arc miitually transversal, 
they do not remain connected. Indeed when, under iteration by Hs, a segment 
of unstable manifold falls over more than one component of the singularity set 
(7, the image is cut into possibly countably many disconnected pieces, as in 
figure 5. We therefore proceed as follows: first, we confirm that exponentially 
lengthening segments remain, even in the presence of cutting; second, we discuss 
images of connected pieces who have grown to the size of Si, for small i; finally, 
we consider Si for large i. 

Observe that the largest eigenvalue of DH is \/5 > 2 and write and 
Ihi') for the vertical and horizontal lengths of a line segment respectively. Then, 
letting 7 be a connected segment of H'g{'j'^{z)), we have that IviHsij)) > 
Hence if i?s(7) lies in exactly two elements Sj, then it contains a connected 
segment 7 such that Ivij) > Sly{'y), for some 6 > 1. Thus Hg{'y{z)) always 
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(a) Initial segment 7 (shown (b) Hg{'y). 

zoomed). 




(c) i?|(7). (d) Stable manifold of c = (1/2, 1/2) 

intersecting Hg{'y). 



Figure 5: Images of a sample initial segment of unstable manifold 7. 5(a) shows 
the initial segment as a black line with a shown in grey, plotted in [0.9, 1] x [0, 0.1] 
for clarity. 5(b) shows the image of 7 under Hs, with the following iterate shown 
in 5(c). In 5(d) we show the hyperbolic fixed point (for Hs) c with its stable 
manifold intersecting a connected component of Hg{{-f)). 
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contains, for each n, an exponentially growing connected segment, at least until 
a segment is cut into three or more pieces. If this occurs then a connected 
segment 7 has spanned one of the elements of Hs from side to side. 

Consider first the region outside J(e), as shown in figure 4(a). and its comi- 
terpart near q. This consists of a finite number of regions S^. It is tedious 
but not difficult, due to the finiteness of the problem, to show that such an 
(Si-spanning segment 7 for each i has the property that a connected segment of 
Hg{j) stretches from y = to j/ = 1 for some finite k. This can be verified by 
considering the images of the endpoints of 7 (these lie on lines of a described 
in appendix A). We show an illustrative example for a typical small initial 7 
in figure 5. We begin in figure 5(a) with a segment assumed to have grown to 
span some 5,. Since 7 lies entirely within an element 5, of cr, Hsij) is again 
connected, and has grown by a factor A'. Careful counting will reveal that in 
this case Hs{{^)) lies over nine elements of a. 5(c) shows the following iterate, 
which consequently contains nine connected components, and which contains a 
connected segment joining y = to y = 1. In this figure that grey backdrop is 
the set a~^, illustrating how connected components of unstable manifolds are 
contained within the singularity set for Hg^. Finally in 5(d) we show the hy- 
perbolic fixed point (for Hs) c with its stable manifold intersecting a connected 
component of if|((7)), and indeed a connected component of all future iterates 

According to [5], the obstacle in demonstrating the Bernoulli property is 
usually the problem that a segment of unstable manifold may be cut into count- 
ably many pieces. In this dynamical system such a phenomenon occurs in 
J(e), or its counterpart near q. Consider now a segment 7 which connects 
the top and bottom edges of Sn, for some sufficiently large n. It has height 
^v{l) > 4(J_i) - 4k = 4n(n-i) - Under a single iterate of Hs, by definition, 7 
undergoes n iterates of F, producing 7' which satisfies ^^(7') = /j,(7) > ^^^n^i^ 
and lh{l) > 2(„-i) ■ ttiis point we have a segment 7' of length ^ ^ which 
now enjoys exponential growth until leaving the component of 6*1 adjacent to 0, 
whereupon it has grown sufficiently that its next iterate produces an S'i-spanning 
set, and the procedure above applies. Considering an identical argument for 
Hg^ this is sufficient to show that there exists to, n € N large enough that 

lH-sr{z)f^Hs"'^\z')^$. 

1 

To demonstrate (4), i.e., that the the equation above holds for a/i sufficiently 
large m and n, consider c = (1/2, 1/2). It is a hyperbolic fixed point for Hs (it 
is a period three point for H). From (4), for all z,z' € S we have 

/iJg7"(^) n 7^(c) ^ m-s'^i'(z') n 7"(c) ^ 0. 

1 The intersections are transversal (see figure 5(d)), so the inclination lemma 
says that successive i?s-images of 7"(z) accumulate on 7"(c), and successive 
ffg^-images of ^^{z) accumulate on 7*(c). Thus (4) follows from (4). 

□ 
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5 Local expansion factors 



In this section we consider the one-step growth condition introduced by Chernov 
and Zhang [8] and related in Section 2. Hs, it transpires, is not itself expansive 
enough and wc arc led to consider Hg instead. 

Let be a global unstable manifold (we suppress the superscript u here 
for case) for let Wi be the connected components of W\(j'^ and let = 
min{A(z) : z G Wi}. Xi is the minimal local expansion factor of Hg on Wi. The 
piece- wise linearity of Hg ensures that A is in fact constant on Wj. Let \W\ 
denote the length ofW. In this section we prove the following: 

Theorem 5.1. We have 




I where the supremum is taken over all such unstable manifolds. 

Let E„ denote those points z e S* for which Hg{z) = The sets 

S„ are to Hg what the sets 5„ arc to Hg. In fact, for large n, 5„ C S„. This 
follows from Lemma 3.2, or by direct construction of cr^. See Figure 6. 

Fix a small £ > 0, let -Be(p), -Be (5) be corresponding neighbourhoods of p 
and q and let Bs{p U q) be their union. We first show that it is sufficient to 
consider W close to either poT q (recall that these are the corners of S at which 
singularities for Hg and its powers accumulate). 

Lemma 5.1. IfWd B^ipUq) = then (5.1) holds. 

Proof. The singularity set cr^ is shown in Figure 6; a^\B^{p U q) consists of a 
finite number of line-segments. Any intersections W f] a'^ 7^ are transversal 
because W and ct^ have all tangent vectors in unstable and stable cones respec- 
tively. Thus a vanishingly short W can cross at most the largest number of 
line segments as meet at any point, which by inspection is two. 

Such a W\cr'^ has at most three connected components Wi so that the sum 
(5.1) contains at most three terms. The maximum value for A,^^ occurs in Si 
and is the largest eigenvalue of {DG-DF f, i.e. (3-^2^2)^ Thus 3/(3-^^2)2 ^ ^ 
is an upper bound for the one-step expansion factor. □ 

It remains to consider W f\ Bs{p q) ^ or, without loss of generality, 

W C B^{pUq). Here singularity line-segments accumulate, with the consequence 
that vanishingly short unstable manifolds W may intersect many of them. 
Our next result describes the expansion factor on S„ for large n. 

Lemma 5.2. If Wn C S„ and A„ is the associated local expansion factor then 
Xn ~ 24n. 

Proof Xn is given by the largest eigenvalue of DHg. It follows from Lemma 3.2 
that for large n there are just four possible forms for DHgi 

DC'-DF-DH, DG-DF'' -DH, DH ■ DG'' ■ DFDH ■ DG ■ DF"", 
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(a) The singularity set 



(b) Sketch of the part of near ■p. 



Figure 6: Illustration of a^, the singularity set for In (a) we show itself, 
to give an idea of the overall structure, without the accumulation of components 
near the corners of the square. Also shown arc the lines L (solid line), V (dashed 
line) and the unstable manifold of p (dotted line). In (b) we show a sketch of tr^ 
in the neighbourhood of p. Here the regions (J and (J accumulate on p. 
The shaded region between the lines L : y = [1 ~ x) /2 and L' : y = 2(1 — x)/5 
effectively contains a copy of a (see appendix A). The light shaded region 
surrounding the dotted line is the component of tr^ which contains the unstable 
manifold of p. 
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where DH = DG ■ DF. In the first case A„ is the largest eigenvalue of 

/ 5 12 
\^ lOn + 2 24n + 5 

given by 12n + 5 + V 144n^ + 120n + 24 24n. The remaining cases are similar. 

□ 

Consider an unstable manifold W C B^(p), having all tangent vectors in C"*". 
A priori Lemma 5.2 poses a problem as follows. By taking W as close to vertical 
as is required and by insisting that W intersects as small a neighbourhood of p 
as is required, W can intersect arbitrarily many sets S„, making the associated 
one-step growth factor J2 l/24n arbitrarily large. This is clearly incompatible 
with (5.1). 

The key to overcoming this seeming difficulty is to observe that W cannot 
be both 'close to vertical' and 'close to p' simultaneously. The following lemma 
formalises this idea. 

Lemma 5.3. Let (u, w)^ be any tangent toWc Bg{p[Jq). Then v/u — )• l + \/2 
as s ^ 0. 

Proof. The result is another consequence of Lemma 3.2. li W C Bg{p U q), 
z £ W r\T,n and e is 'small' then either n = 1 or n is 'large'. 

First suppose n is large, then by Lemma 3.2 at least rn ~ win n, k > 0, of 
the immediate iJ^^-images of z are in Ei. Let {u',v')'^ G be a tangent to 
Hg^'^iW) at Hg^'^iz) = The tangent to W at z is then given by 

{u,v)'^ = {DG ■ DF)'^"'{u',v')'^. As £ 0, so m oo, and the tangent to W 
approaches the imstablc eigenvector of DG ■ DF, i.e. v/u ^ 1 + \/2- 

Conversely suppose z G fl Si. Although Lemma 3.2 doesn't apply imme- 
diately, the proximity of to p or g has the same consequence that some number 
of the immediate i/^ ^-images of z are in Si, and from here the argument is as 
above. □ 

To establish (5.1) it is enough (Lemma 5.1) to consider W approaching p (or 
q, but without loss of generality we focus on the former). Such a W can intersect 
arbitrarily many sets S„ but our control over its orientation (Lemma 5.3) and 
knowledge of the growth factor on S„ (Lemma 5.2) will lead to a bound on the 
corresponding sum J2 ^n^- The next lemma establishes such a bound. 

Consider the conical region of B^{p) bounded by L (satisfying y = (1 — x)/2) 
and hy X = 1. Let S^ be the component of S„ in this region. We consider the 
connected component of W confined to this region; the case of more general W, 
establishing the theorem, follows the lemma. 

Lemma 5.4. Let W be a connected component of an unstable manifold, of 
length 5 > 0, with ends on L and x = 1. Then 

lim inf V = In f 3 + 2V2) . 

n 
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Proof. Let W O L — (1 — 2yo,2/o) for some small yo > 0. Lemma 5.3 says that 
W intersects a; = 1 at y = t/i « (3 + 2-\/2)yo, with equality in the limit 5 — >■ 0. 

The 'lower' boundary of S^, which is also the upper boundary of ^'n+i^ 
intersects x = 1 a,t y = l/2n and intersects L at y = l/2(n — 1). Let N be the 
unique integer so that 

1 1 

2N "^^^^ 2{N-1)' 

W intersects S'^ but not S'^ for any n < N. The asymptotic relationship 
between j/i and yo says that there is an integer M depending on N so that W 
intersects E^, W does not intersect E^ for any m > M, and M (3 + 2\/2)A''. 
The limit N ^ oo corresponds to 5 — ^ 0, thus 

M{N) r(3+2V2)JV] 

liminfyA^i= lim V A'^ = lim V — = — ln(3 + 2^), 
where we have used the fact that , ^ — In — >■ const. □ 

/-^n—l n 

We now prove the main result of the section. 

Proof of Theorem 5.1. In light of the comments made following Lemma 5.3 we 
let W have end-points on y = and x = 1 and gradient 1 + \/2. Such a W 
intersects four distinct regions where singularities accumulate, and the region 
El. We determine the contribution to liminf^^o suppj^.|^|<^ of each. 

Lemma 5.4 deals with the region IJ E^, bounded by L and x = 1. Now let 
E^ be the connected component of E„ in the region bounded by L' satisfying 
y = 2(1 — x)/b and by y = 0. Let r denote reflection through x + y = 1. Notice 
that 

roF(E;;) = E; 

for each n. Moreover the gradient of is (r o F)-invariant. Thus Lemma 5.4 
applies to r o F (IJ EJ() and, because r o F is invertible, to IJ E^ itself. 

The remaining two regions where singularities accumulate can be dealt with 
in the same manner. For the region adjacent to |J E'^ the appropriate invertible 
transformation \s G o F and for the region adjacent to IJ E" it is G o F o r o f . 

Finally, W crosses Ei where, by Lemma 5.1, \{z) = (3+2\/2)^. We conclude 
that 

liminf sup V Ar^ = ^ + - ln(3 + 2^2) < 1. 

5^0 W:\W\<5'^ ' (3 + 2x/2)2 6 



□ 



6 Proof of the main result 

The following theorem shows that the analogue of (2) is satisfied. 
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Theorem 6.1. There is a set A C S of positive ns-measure and having hyper- 
bolic product structure, and 9 e (0, 1) so that 

^i{zcS: Rtn* (z;iJs, A) > n} = 0(6*"). 

Proof. As described in Section 2 it is not necessary to explicitly construct A. 
Rather we show that Hs satisfies the conditions, essentially due to Chernov [5], 
that were listed. The necessary work was complcitcd in Sections 3, 4 and 5. 
We remind the reader where each result may be found; italics correspond to 
subsection headings of Section 2. 

The smoothness condition concerns the set a defined in Section 3. It is ev- 
idently closed and a countable union of zero-measure line segments. Much of 
the hyperbolicity condition is demonstrated in Lemma 3.1 with the remainder, 
concerning Lyapunov exponents and local invariant manifolds, in Lemmas 4.1 
through 4.3. Existence of an invariant SRB measure is given by Proposition 3.1. 
That it is mixing follows from the stronger Bernoulli property, proved in The- 
orem 4.1. The conditions on distortion bounds and bounded curvature follow 
immediately from the piecewise linearity of H and thus of Hg. Indeed, the ex- 
pansion factor must be constant on local invariant manifolds, which themselves 
are zero-curvature line-segments. Absolute continuity of the foliation follows 
from the result of [13] and the Lemmas 4.1 through 4.3. The condition on the 
structure of the singularity set holds because any unstable curve W intersects 
a transversally and at most countably many times; the only possible accumu- 
lation points are p and q and the rate of convergence along such a sequence of 
intersections is of order 1/n, owing to the structure of a. The one-step growth 
condition was the subject of Theorem 5.1. Finally, /i and differ by a constant 
factor on S so that the statement holds as given, i.e. with /i rather than iig. □ 

Theorem 6.1 shows that Hs : S ^ S has exponential decay of correlations 
for a-Holder observables and constitutes the majority of our work in proving 
Theorem 1.1. 

To conclude the proof we establish (2.1), i.e. with A as above and with 
An = {z & R: Rtn*(^;iJ,A) > n}, 

we show that 

Theorem 6.2. /z(A„) = 0{l/n). 

As mentioned in Section 2 we follow a procedure introduced in [15]. It 
involves treating separately a certain set of 'infrequently returning' points, to 
be defined. For n e N and ^; e i? let 

n 

r{z;n,S) = ^Xs {H\z)) 

be the number of the first n images of z that are in S. Let 6 > be a constant 
(to be fixed shortly) and define 

Bn,b = {z & R: r{z; n, S) > 6 Inn}. 
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Bn,b contains those points returning to S at least 6 In n times within n iterations. 

Lemma 6.1. /i(A„ n Bn,b) = 0{l/n). 

The proof is due to [15, 8] but is included for the sake of completeness. 

Proof. Let z G And Bn^ and let i = Rtn(2:; H, S). Clearly < i < n (in fact 
i < n — b\nn, but the weaker bound will suffice). Prom the definitions of An 
and Bn,b we have 

Rtn* {H\zy,Hs, A) >b]iin. 

and so 

n-l 

A„nB„,bC U H-'{z' G S -.Rtn* {z';Hs, A) > bin n}. 

Theorem 6.1 now gives 

niAnDBnfi) < nn{z' € S : Rtn{z';Hs,A) > blnn} = O (n6l^'"") . 
Taking b > -2/ In 61 > gives 6*""" < n''^ and thus the result. □ 
Let us now fix 6 > —2/\n6 > as above. For z G R and n e N let 

iVmax(z, ") = max {Rtn {H\z);H, S) : Q i n] . 

^raa.yi{z,n) is the largest interval either until we first enter S, or between con- 
secutive returns to S. The main step in our proof that /x(A„\i?„^f,) = 0{l/n) 
is to show that, away from Bn,b, -^max grows linearly. 

Lemma 6.2. There is a constant p > such that if r{z;n, S) ^ blnn then 
Njm,^{z,n) > /3n. 

Proof. Suppose that r{z:n, S) ^ blnn. N-a-iayi{z,n) is minimised when the total 
return-time n is distributed as evenly as possible between the 6 Inn returns. For 
small n the mean return-time n/blnn gives a lower bound. 
Now suppose that 



6 Inn 

where fc, K arc the constants of Lemma 3.2. Some returns necessarily land in 
Sn, N ^ Ke'^, and so adjacent returns are to ^i. Thus evenly distributing the 
total return-time n between the 6 Inn returns is inconsistent with Lemma 3.2. 
In this case Njjia,^{z,n) is minimised by an itinerary of the form 

5i, Sn^,Si, S'jVa, i^i, 5jv3, ... 

where each Ni is approximately 2n/b\nn— 1 < A^max(-z, n). 
If additionally 

^-l^Ke'^ 
bmn 
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then the above arrangement into pairs is also inconsistent with Lemma 3.2. Here 
■^max is minimised by an itinerary 

■ ■■,Si, Si,Sni, Si, Si, Sn2iSi, Si,Sn3, ■■■ 

where each Ni is approximately 3n/61nn — 2 ^ ■^max('2) n). 
In general if j G N and 

IjP- - {j - 1) > Ke^'' 
In n 

1 then 

bmn 

1 For n e N let J{n) be the largest integer j for which (6) holds, then 

i^^±^ - J < ife(^+^)^ 
olnn 

It follows that 

l-^ < fi^e(^+i)'= + j)< Ke^'+^^K 

bmn J+IV / 

1 The first inequality in (6) is a rearrangement of the previous displayed equation. 
The second follows easily from the assumption that Ke^"^"'"^^'^ > 1, which holds 
for all sufficiently large J, i.e. for all sufficiently large n. 

Taking logarithms on each side of (6) and rearranging gives 

J + 1 > ^ (In n — In In n — k) 

k 

where k = IniiT + lnfe is constant. By assumption (6) holds with j = J therefore 
(6) gives 

,^ , , n / Inlnn k \ 1 , , , ^ 

A^.ax(.,n) >^(l-^-i^J-^ (Inn- Inlnn- + 1. 

Clearly Ny^ax/n — )-l/A;6>0asn— J- oo, at a rate that is independent of z. 
Hence the lemma holds for some < ^ < 1/kb. □ 

Lemma 6.3. /i(A„\B„,b) = 0{l/n). 

Proof. We prove the sufficient result ^{R\Bn^b) = 0{l/n). By Lemma 6.2 

fi{R\Bn,b) = iJ,{z G R : r{z;n, S) ^ blnn} 
< IJ,{z G R: N^a;,r{z,n)^ I3n} . 

Recall that Nynax{z,n) is either the largest N such that z lands in Sn, 
or the number of iterations taken to first enter S. If it is the former then 
W{z) e {Sn ■ N > /3n} for some < i < n. We have seen, in the proof of 
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Lemma 4.2, that /i(S'„) = 0{l/n^), and thus /x{5jv : |A''| > M = ^'(l/?^^), 
and thus 

M (ii\B„,6) < M /^-'{■^jv : |A^| > ^n} j = 0(l/n). 

Conversely, notice that z ^ S and let H^{z) E S he the smallest such j. 
For n G N let S'-n be the elements of the partition of S induced by Hg^ 
(these arc analogous to the sets Sn induced by Hg). Then H^{z) e S-m and 
M > Afmax(2, n) > /3n, and the argument proceeds as above. □ 

Theorem 6.2 follows immediately from Lemmas 6.1 and Lemmas 6.3. This 
completes the proof of Theorem 1.1. 

Finally we remark on the changes required to our proof in order to accom- 
modate different annuli P and Q. The differences most obviously manifest 
themselves in the values of certain constants. Lyapunov exponents in particular 
will vary, impacting the various constants in Sections 3 and 4. The exact val- 
ues however are unimportant to the arguments, and important properties, e.g. 
positivity, will not change. Geometric features such as the precise structure of 
a and the relative sizes of the sets 5„ will also change, but again the important 
features, such as the general structure of a and the boimd /x(S'„) = 0(l/n^), 
remain. Lastly, if local growth rates are weaker than in the map considered, we 
might need to consider a higher iterate of Hs in Section 5, but this introduces 
no further difficulties than have presently been overcome. 



7 Correlation decay in other linked-twist maps 

We have shown that the rate of decay of correlations for a large class of linked 
twist maps is polynomial. However, note that if both of the annuli are thickened 
until they are equal to the entire torus, the map H becomes the hyperbolic 
toral automorphism known as the Arnold Cat Map, which is well-known to 
be exponentially mixing. The transition from non-uniformly hyperbolic linked- 
twist map to uniformly hyperbolic Cat Map deserves further study. 

Similarly, increasing the wrapping number of the twists (that is, taking at 
least one of j, k > 1) only serves to enhance the mixing, yet the behaviour at 
the boundary remains linear. In this case we expect only minor modifications 
to produce identical results. An interesting case arises if the twists / and g are 
permitted to be nonlinear, yet still monotonic. As in [3, 20], H is still Bernoulli, 
although the behaviour at the boundaries may now be different. Again, we 
expect the dominant behaviour to be sub-exponential in this situation. 

If exactly one of j and k is allowed to be a negative integer, the situation is far 
more complicated. Although in this case H can still be shown to be Bernoulli, for 
certain choices of / and g, the proof relies on an intricate geometrical argument 
due to [20] . The question of its rate of correlation decay is still open. Likewise 
the LTMs defined on planar annuli of [26] are Bernoulli [22], but the methods 
in this article would need significant adaption. 
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(a) 



(b) sf^ 



(c) F- 1 {Sf ^ ) (d) n F- 1 {S[^^ ) 



Figure 7: The construction of the set ^i, consisting of all points of S which 
return to 5 under one iterate of H. Sets Si for larger i can be considered in the 
same way, and their union forms the partition in figure 2(c). 
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A Structure of the singularity sets 

In this appendix we give details of the construction of the singularity sets for 
Fg, Gs, Hs and Each set consists of the pre-iniages of the boundary of S 
under the return map in question, and which partition S into distinct regions 
which take different numbers of iterates to return to S under the map {not the 
return map) in question. 

For example, the lower half of the singularity set for Fs consists of a sets of 
lines connecting 2y ^ 1 — x with x = 1. Each line in this set meets 2y = 1 — x 
at ((n — 2)/(n — 1), l/2(n — 1)) and meets x = 1 at (1, l/2n), for each n > 2. 
These lines are just the lines which take n iterates of F to be mapped into the 
line a; = 0, and hence to return to S. Similarly, the upper half of this singularity 
set contains a set of lines connecting 2y = 2 — x with x = 0, which are the lines 
which takes n iterates to be mapped into x = 1. These sets of lines accummulate 
on q and p respectively. See figure 2(a). 

The singularity set of Gs, shown in figure 2(b), is exactly analogous. This 
set is just the singularity set for Fs reflected about the line y = 1 — x. Again, 
constituent lines accummulate on p and q. 

The singularity set for Hs, denoted a, can be constructed in a similar way. 
We are also interested in the sets Sn = {z € S : 'Rtn{z; H, S) = n} so, analo- 
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gously, we define the sets 

5^^) = {zeS: Rtn{z; F, S) = n} 
S'^c) ^ {zgS: Rtn{z; G, S) = n} 

These sets can be easily discerned from figures 2(a) and 2(b), with and 

(G) 

Sn accumulating on p and g as n — > oo. Now by definition of Hs, the sets 5„ 
have the property that a point 2; € 5„ if and only if 2: € Sj^^ nF~^sl^^), where 
j + k — n + 1. Wc illustrate this statement by constructing explicitly, as 
shown in figure 7. Points in the shaded region of figure 7(a) return to S under 
j = 1 iterates of F. Points in the shaded region of figure 7(b) return to S under 
A; = 1 iterates of G. Thus to return to S under a single iterate of H, a point 
must lie in both , and the pre- image under F of Si'^ (shown in figure 7(c)). 
This intersection is shaded in figure 7(d). 

A similar statement can be made for a^. A point z e I]„ if and only if 

z e S^.p n i^-^'^(4?) n F-^'HG-'=i(s]f )) n f-^HG-'^(f-^'m4?)))' (3) 

where ji + fci + j2 + fc2 = n + 3. Thus the statement that the shaded region of 
figure 6 effectively contains a copy of a can be understood by the fact that this 
shaded region is part of Si, and equation (3) must have ji — ki = 1, leaving j2 
and k2 to satisfy j2 + k2 = n + 1, just as in the construction of a. 
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